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ABSTRACT 

Measuring the 3D distribution of mass on galaxy cluster scales is a crucial test of the ACDM 
model, providing constraints on the behaviour of dark matter. Recent work investigating mass 
distributions of individual galaxy clusters (e.g. Abell 1689) using weak and strong gravi- 
tational lensing has revealed potential inconsistencies between the predictions of structure 
formation models relating halo mass to concentration and those relationships as measured in 
massive clusters. However, such analyses employ simple spherical halo models while a grow- 
ing body of work indicates that triaxial 3D halo structure is both common and important in 
parameter estimates. Though lensing is sensitive only to 2D projected structure and is thus 
incapable of independently constraining 3D models, the very strong assumptions about the 
symmetry of the lensing halo implied with circular or perturbative elliptical NFW models 
are not physically motivated and lead to incorrect parameter estimates with significantly un- 
derestimated error bars. We here introduce a Markov Chain Monte Carlo (MCMC) method 
to fit fully triaxial models to weak lensing data that gives parameter and error estimates that 
fully incorporate the true uncertainty present in nature. Using weak lensing data alone, the 
fits are sensitive to the Bayesian priors on axis ratio; we explore the impact of various general 
and physically motivated priors, and emphasise the need for future work combining lensing 
data with other data types to fully constrain the 3D structure of galaxy clusters. Applying the 
MCMC triaxial fitting method to a population of NFW triaxial lenses drawn from the shape 
distribution of structure formation simulations, we find that including triaxiality cannot ex- 
plain a population of massive, highly concentrated clusters within the framework of ACDM, 
but easily explains rare cases of apparently massive, highly concentrated, very efficient lensing 
clusters. Our MCMC triaxial NFW fitting method is easily expandable to include constraints 
from additional data types, and its application returns model parameters and errors that more 
accurately capture the true (and limited) extent of our knowledge of the structure of galaxy 
cluster lenses. 

Key words: gravitational lensing - cosmology: theory - dark matter - galaxiesxlusters: gen- 
eral. 



1 INTRODUCTION 

Galaxy clusters are ideal laboratories in which to study dark mat- 
ter, being the most massive bound structures in the universe and 
dominated by their dark matter component (~ 90%). Constrain- 
ing the clustering properties of dark matter is crucial for refining 
structure formation models that predict bo th the shapes of dar k 
matter halos and their mass function (e.g. iNavarro et al.l ( Il997h : 
[fiahcall et al. ( 2003); Dahle (2006)). Several methods are used to 
measure galaxy cluster dark matter profile shapes and halo masses 
on a range of scales, including X-ray and Sunyaev-Zeldovich (SZ) 
studies, dynamical analyses, and gravitational lensing. However, 
all of these methods require simplifying assumptions to be made 
regarding the shape and/or dynamical state of the cluster in or- 
der to derive meaningful constraints from available data. Crucially, 
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while most parametric methods typically assume spherical symme- 
try of the halo, observed galaxy clusters often exhibit significant 
projected el lipticity and halo s in CDM structure formation simu- 
lations (e.g . iBett et ai] ( l2007l) (using the Millennium simulation); 
IShaw et all ( l2006h ) show significant triaxiality in cluster-scale ha- 
los, with axis ratios between minor and major axes as small as 
0.4. Understanding and accurately incorporating the impact of this 
physical reality on cluster mass and parameter estimates is crucial 
for accurate comparisons between measured cluster properties and 
model predictions. 

In addition to the determination of masses, most cluster pro- 
file fits are carried out in the hope of either supporting or refut- 
ing the universality of the NFW profile and thus testing the CDM 
paradigm. The NFW profile is typically parameterised by an ap- 
proximate virial mass i\f200 and a concentration parameter, C, and 
simulations predict a strong correlation between the two. For a clus- 
ter of M — 10^^ Mq, C ^ 4. However, several authors (e.g. 
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iLimousin etal] ( l2007h : [Xneib et alj l2003h : lOavazzi et alj Jiooj) ') 
have recently reported gravitational lensing results in the very low 
probability tail of the predicted distribution; notab ly, in a combined 
weak and strong lensing analysis of Abell 1689. lBroadhurst et alJ 
( l2005h report a concentration parameter of C ~ 10, when C ~ 4 
is expected. While a few such results are not damning, especially 
given the very complex, likel y not relaxed , struc tures of the lens- 
ing halos (see e.g. the work of Lok as et al.l ( |200^ on A1689), it is 
nonetheless of interest to investigate how possible future discrep- 
ancies between observations and the predictions of ACDM should 
be interpreted. 

Efforts to understand the impacts of triaxiality in gravitational 
lensing and its potentia l role in exp l aining apparent discrepancies 
with CDM began when lOguri et al.l ( 1200^) applied a fully triaxial 
NFW model to the shear map of Abell 1689 to find that it is con- 
sistent with 6% of cluster-scale halos. These continued as lGavazzil 
showed that a triaxial NFW can reconcile parameter values 
derived from observations of the cluster MS2 137-23 to predictions 
from N-body simulat ions. 

More generally, ICorless & KiriS ( l2007h demonstrated in the 
weak lensing regime that neglecting halo triaxiality in parame- 
terised fits of NFW models to NFW haloes with axis ratios signif- 
icantly less than one can lead to over- and underestimates of up to 
50% and a factor of 2 in halo mass and concentration, respectively. 
While extreme cases of triaxiality are rare, such haloes can be much 
more efficient lenses than their more spherical counterparts, espe- 
cially when in configurations that hide most of the triaxial shape 
along the line of sight. Further, even haloes with less extreme axis 
ratios are inaccurately fit by spherical models. We expect the vast 
majority of galaxy cluster scale dark matter haloes, and so it is im- 
portant to include this expectation in model fitting. 

In the past, triaxial models have not generally been fit to lens- 
ing data because they cannot be well-constrained. A simple NFW 
triaxial halo model has six free parameters (concentration, mass, 
two axis ratios, and two orientation angles) as opposed to the two 
parameters of the spherical model (mass and concentration only), 
and even combined weak and strong lensing data is not enough to 
meaningfully constrain so many parameters. This is sometimes due 
to the limited depth of currently available observational data, but 
more importantly to the intrinsic limitations of lensing. Because 
lensing is affected only by the the projected surface density and 
shear of the underlying mass distribution, it is inherently impossi- 
ble to fully constrain a three-dimensional structure without impos- 
ing strong priors on the shape of the halo or supplementing lensing 
data with other data types more sensitive to line-of-sight halo struc- 
tures, such as dynamical studies. 

Despite these difficulties, since triaxiality has been convinc- 
ingly demonstrated to be an important factor in model parameter es- 
timation in lensing analyses, it is necessary to find a way to directly 
include it in NFW fits to galaxy clusters. For individual clusters it 
is crucial, if claims regarding the validity of the ACDM paradigm 
based on NFW parameter estimates are to be meaningfully eval- 
uated. It is also of importance across populations to obtain more 
accurate distributions of galaxy cluster parameters, and especially 
in measuring the galaxy cluster mass function - if the mass func- 
tion is signficantly sloped as expected, even the best-case scenario 
of a symmetric scatter of mass estimates due to neglected triaxiality 
would lead to an asymmetric shift in the calculated mass function, 
as there are more low mass halos to shift up in mass than there are 
high mass halos to shift down. 

To that end, we present here a Bayesian Markov Chain Monte 
Carlo method to fit fully triaxial NFW halos to weak lensing data. 



This method flexibly combines weak lensing data with prior prob- 
ability functions on the model halo parameters to return parameter 
and error estimates that reflect the true uncertainties of the problem. 
Importantly, it also allows for the straightforward and statistically- 
robust inclusion of additional constraints from other data sources 
such as observations from SZ, X-ray and spectroscopic surveys. 



2 LENSING BY TRIAXIAL NFW HALOS 
2.1 Weak Lensing Background 

Weak lensing distorts the shapes and number densities of back- 
ground galaxies. The shape and orientation of a background galaxy 
can be described by a complex ellipticity e", with modulus je°| = 
(1 — b/a)/{l + b/a), where b/a is the minormajor axis ratio, 
and a phase that is twice the position angle 4>, e" = je^je^*'^. 
The galaxy's shape is distorted by the weak lensing reduced shear, 
g = 7/(1 — k), where 7 is the lensing shear and k the convergence, 
such that the ellipticity of the lensed galaxy e becomes 



+ 7 



(1) 



1 = 

in the limit of weak deflections. The distributions of ellipticities for 
the lensed and unlensed populations are related by 



(2) 



assuming a zero-mean unlensed population, the expectation values 
for the lensed ellipticity on a piece of sky is < e >= g « 7. 
This is the basis for weak lensing analysis in which the shapes of 
images are measured to estimate the shear profile generated by an 
astronomical lens. 

Lensing also changes the number counts of galaxies on the 
sky via competing effects; some faint sources in highly magnified 
regions are made brighter and pushed above the flux limit of the 
observation, but those same regions are stretched by the lensing 
across a larger patch of sky and so the number density of sources is 
reduced. Thus the number of sources in the lensed sky n is related 
to that in the unlensed background no and the slope of the number 
counts of sources at a given flux limit a by n = noj-L"~^, where 
fi is the lensing magnification fj,~^ = (1 — k)^ — |7|^. A full 
description of these effects is given in lCanizMe 

We now describe the characteristic behaviour of the conver- 
gence n and shear 7 of triaxial NFW halos. 



2.2 Triaxial NFW 

A full parameteri sation for a triaxial NFW halo is given by 
Ijing & Sutd ( I2OO2I) (herein JS02). They generalise the spherical 
NFW profile to obtain a density profile 



p(i?) = 



ScPc{z) 



R/Rsil + R/RsV 



(3) 



where 5c is the characteristic overdensity of the halo, pc the critical 
density of the Universe at the redshift z of the cluster, Rs a scale 
radius, R a triaxial radius 

Y'S y2 

fl' = — + ^ + — , (a^b«:c=l), (4) 
b^ 

and a/c and b/c the minormajor and intermediate:major axis ra- 
tios, respectively. In a different choice from JS02 we define a triax- 
ial virial radius R200 such that the mean density contained within 
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an ellipsoid of semi-major axis -R200 is 200pc such that the con- 
centration is 



C 



-R200 



(5) 



the characteristic overdensity is 
200 



5c = 



log(l + C)-^' 



(6) 



the same as for a spherical NFW profile, and the virial mass is 
SOOtt 

M200 = abH2ooPc- 



(7) 



In past lensing studies, an effective spherical virial mass and 
concentration have often been employed instead of these triaxial 
definitions, in part to allow more straightforward comparison to 
populations of structures extracted from simulations using spher- 
ical halo boundaries. The effective spherical virial radius is defined 
as the radius r2oo at which the mean density within a sphere of that 
radius is 200 times the critical density, and the effective spherical 
virial mass the mass within that sphere: 



m,2oo = (8007r/3)r2ooPc 



(8) 



We further define the effective spherical concentration Caph of a 
triaxial halo as the ratio of the effective spherical virial radius to 
the geometric mean of the triaxial scale radii rs = Ra (abc)^^^: 



sph 



■ r 200 /rs 



(9) 



We test the performance of our fitting method in recovering these 
effective values for the purposes of comparison with previous work; 
however, we stress that a fully triaxial definition of virial mass bet- 
ter represents the true shape of dark matter halos, and is well mo- 
tivated by ellipsoidal collapse models that predict asymmetric col- 
lapse to stop at the same enclo sed density as does spherical collapse 
dSheth, Mo, & TormeiJ JioOlh ). We therefore prefer the fully triax- 
ial parameterisation in this and future work as the more physically 
motivated choice. In practice, the effective spherical mass and con- 
centration are always less than and vary only moderately from their 
fully triaxial counterparts, by up to 10% for very triaxial halos and 
most often far less. 

The triaxial halo is oriented with respect to the observer's line 
of sight by angles 9 and 0; randomly oriented halos are distributed 
uniformly in <j) and sin^. A more deta iled description of this p aram- 
eterisation and its benefits is given in lCorless & KingI ( l2007t) . 

The f ull derivation of the l ensing properties of a triaxial halo 
IS given bv'Oguri, Lee, &Sutol ilOQJi (herein OLS), and we sum- 
marise and extend some of that work in Appendix [A] 



2.3 Weak Lensing Simulations 

The main body of simulations is carried out for a field 7.5' in ra- 
dius, with a background source density no = 30/arcminute^ (Pois- 
son noise is accounted for), typical of ground-based observations. 
Although wide-field imaging of clusters with fields of half a degree 
is routinely possible, there are arguments that errors due to large- 
scale structure along the line-of-sight become more important as 
the shear due to the cluster itself diminishes with distance from 
the centre (e.g. Hoekstra (2003)). In any case, the general trends 
will hold for larger fields. A catalogue of randomly positioned and 
oriented galaxies with intrinsic shapes e" drawn from a Gaussian 
distribution with dispersion = a^i + = 0.2 in the mod- 
ulus je"! is placed at redshift z = 1. This catalogue of background 




0.6 
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Figure 1. The prior probability functions for axis ratios b (solid lines) and 
a/b (dashed lines). The bold lines show the Shaw prior, the o-dotted lines 
the Axis Gaussian prior generaUsed from Bett et al., and the Ughtest Unes 
the Flat prior. 




Figure 2. The top {lower} panels show isoconvergence contours for a 
prolate {oblate} halo of M200 = lO'^^ Mq, C = 4, with axis ratios 
a = b = 0.4 {a = 0.4, b = 1.0}; the left-hand panel shows the halo ori- 
ented with the odd axis along the line of sight (LoS), the right-hand panel 
shows the halo with odd axis in the plane of the sky (Plane). The thick solid 
lines show the limits of the aperture from which weak lensing data is taken, 
and the dashed line shows the Rs ellipse for each projection (note that Rs 
is scaled by the axis ratio in each direction, so that when looking at the mi- 
nor axis of a triaxial halo, the apparent scale radius is a times the Rs value 
of the halo). The lowest contour corresponds to k = 0.02 and subsequent 
contours each increase in k by a factor of 2. 
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Figure 3. M200 — C confidence contours at 68% and 95% for the triaxial NFW fit to four highly triaxial lenses under various priors. These ai'e the four 
extreme cases noted in our first triaxiality paper: symmetric oblate and prolate halos with Af200 = 10^^ M0 and C = 4 and axis ratios {a = 0.4, b = 1.0} 
and {a = 6 = 0.4}, respectively, aligned in Line-of-Sight and Plane of the Sky orientations. The black stars show the true parameters of the underlying lens. 



galaxies is lensed through a model lens of choice placed at red- 
shift 2 = 0.18 (the redshift of Abell 1689), at which the width 
of the field is ~ 1900 kpc/h. Thus our choice to place all sources 
on a sheet at z = 1 is justified by the low redshift of our fidu- 
cial lens; only for higher redshift lenses that are in the heart of the 
reds hift distribution is the di stribution of source redshifts impor- 
tant ('Sei tz & Schneideilll997h . The background galaxies are lensed 
according to Equation [l] and the number counts are reduced as pre- 
scribed in section im taking the slope of the source number counts 
in flux to be A\ogN / A \o%S = a = 0.5 (correspond i ng to a slope of 
0.2 in magnitude as in iFort. Mellier. & Dantel-For^h997l) ). Galax- 
ies located within 1' of the cluster centre are removed from the 
analysis to avoid the strong lensing regime at the centre of the 
cluster (in any case background galaxies near the cluster centre 
would be mostly obscured by cluster members in observations). 
Throughout we assume a concordance cosmology with flm =0.3, 
Ho = 70 km s^^ Mpc~^, and a cosmological constant Qa ~ 0.7, 
and a typical massive cluster of triaxial M200 ~ 10^^ Mq and 
C = 4. At the redshift of our lens z = 0.18, 1 arcminute corre- 
sponds to ~ 127 kpc//i. 



3 MARKOV CHAIN MONTE CARLO FITTING 

Markov Chain Monte Carlo (MCMC) methods are hugely valuable 
in under-constrained or highly-degenerate fitting problems; they are 
employed widely in CMB analyses and are rapidly growing in pop- 
ularity in lensing work (see e.g. [Marshall (2006), Julio et al. (2007]), 
ICorless. Dobke. & Kind ( I2OO8I) ). This family of methods employ a 
"guided" random walk that returns a sample of points representa- 
tive of the posterior probability distribution; the probability of a 
certain region of parameter space containing the true model is di- 
rectly proportional to the density of points sampled in that region. 
From the distribution of sample points the full posterior probability 
distribution is obtained, which is easily and directly marginalized 
over to obtain fully marginalized mean most-probable parameter 
estimates for all parameters. These most-probable parameter values 
reflect the full shape of the posterior probability distribution with- 
out assumption about the shape of the error distribution, and taking 
full account of the prior. Such methods have exploded in popular- 
ity recently, and there ar e several excellent references describin g 
the method in detail, e.g. lLewis & Bridld ( l20o3) . lMacKavl ( l2003h . 
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Figure 4. ID parameter distributions for the triaxial NFW fitted to weak lensing by four highly triaxial halos under various priors. These are for the four 
extreme cases noted in our first triaxiality paper: symmetric oblate and prolate halos with M200 = 10^^ Mq and C = 4 and axis ratios {a = 0.4, b = 1.0} 
and {a = b = 0.4}, respectively, aligned in Line-of-Sight and Plane of the Sky orientations. The thick black lines show the true parameters of the underlying 
lens. 



Simply put, the sampler functions by stepping through parameter 
space by taking random steps governed by a transfer function, usu- 
ally a simple n-D Gaussian, where n is the number of parameters 
of the fitted model. If the randomly-chosen next step is to a point 
of higher probability than the current position, the step is taken. 
If the next step is to a point of lower probability, the step is taken 
with probability p(new)/p(currcnt). Thus, the MCMC sampler 
spends most of its time in high probability regions, but can move 
"downhill" to regions of lower probability in order to explore the 
entire space and sample all regions of high probability. Crucially, 
this method is able to return a true representation of the full poste- 
rior probability distribution, regardless of the complexities of that 
probability distribution, e.g., tight degeneracies, multiple islands of 
high probability, or a very flat distribution due to poor constraints 
from the data. 



give complete information about the 3-D triaxial structure because 
all lensing behaviour is determined only by the 2-D projected mass 
density and potential. With such weak constraints, a maximum like- 
lihood approach is both impractical and of very limited scientific 
value; the posterior distribution will be very flat with significant de- 
generacies giving poorly constrained maximum likelihood values. 
By exploring the full posterior probability distribution, this MCMC 
method allows the derivation of parameter estimates (and their ac- 
companying errors) that account for the true uncertainties when fit- 
ting parametric models to lensing data. 

To implement the method, we must define the posterior prob- 
ability function, defined in Bayesian statistics as 



p{e\Tv)p{Tv) 
p{d) 



(10) 



Fitting triaxial models with lensing data is an intrinsically 
under-constrained problem, as lensing on cluster scales can never 



where p{6\tv) is the likelihood £ of the data given the model pa- 
rameters (the standard likelihood), p{n) is the prior probability 
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Figure 5. The top panels show the m2oo — Csph confidence contours at 68% and 95% for the effective spherical parameterisation of the triaxial NFW fit to 
the Line-of-Sight oriented highly triaxial prolate and oblate lenses of Figure[3]under various priors. The black stars show the true effective spherical parameters 
of the underlying lens. The bottom panels plot the ID effective spherical parameter distributions for the same lenses and priors, for comparison withE] The 
thick black lines show the true effective spherical parameters of the underlying lens. 



distribution for the model parameters (e.g. a distribution of axis 
ratios drawn from simulations), and p{6) is a normalising factor 
called the evidence, of great value in comparing models of dif- 
ferent classes and parameter types, but expensive to calculate and 
unnecessary for the accurate exploration of the posterior distribu- 
tion. The assignment of priors is often a controversial exercise; we 
will return to this question since priors acquire an increased signif- 
icance in under-constrained problems such as ours. We define the 
log-likelihood function in the standard manner for weak lensing 
following . Schneider. King. & Erben C2000.) and , King & Schneider. 
( l200lh 

^ = -ln£ ^ ^^\np,{e,\g{0,-n)). (11) 

i = l 

where the reduced shear g is calculated using the triaxial conver- 
gence and shear of Eouations lA12l and 11 is a six-element vector of 



the parameters defining the model: triaxial virial mass M200, con- 
centration C, minor axis ratio a, intermediate axis ratio 6, and two 
orientation angles 9 and (j). 

In our MCMC sampler we employ a 6D two-sided Gaussian 
transfer function, and use the covariance matrix of an early run to 
sample in an optimised basis aligned with the degeneracies of the 
posterior. We tune the step sizes of the sampler to achieve an av- 
erage acceptance rate of 1/3 in each basis direction, run three in- 
dependent MCMC chains, started at randomly chosen positions in 
parameter space, for each lens system, and sample the distribution 
space until the standard var(chain mean)/mean(chain var) indicator 
is less than 0.2, indicating chain convergen ce. We util ise the Get- 
Dist package from the standard CosmoMC jLewis fe Bridle 20o3) 
distribution to calculate convergence statistics, parameter contours, 
and marginalized parameter estimates. We employ 40 bins in the 
Gaussian smoothing of the contours, chosen because, for the spher- 
ical case, this smoothing length best matches the contours obtained 
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Figure 6. The full posterior probability distribution for a triaxial NFW fitted to weak lensing by an extreme prolate halo a=b=0.4 in a Line of Sight orientation, 
under a Flat prior. The black contours give the 68% and 95% confidence contours of the marginalized posterior probability distribution, and the shading shows 
the marginalized likelihood values, smoothed over 20 bins (rather than 40 as in other figures) for increased clarity. The solid (dotted) lines in the unshaded 
distributions at the end of each row give the 1-parameter marginalized probability (likelihood) distributions. 



using this MCMC method to those of a standard approach. We 
note that while the angles are defined over a range < 6* < tt and 
Q < 4> < 2n, because of the elliptical nature of the projected den- 
sity contours, they give rise to unique lensing profiles only over the 
range of < 6* < 7r/2 and < < tt. 



3.1 Choice of Priors 

As noted above, the choice of priors used on the parameters of the 
triaxial NFW model is very important because of the inherently 
under-constrained nature of the problem. The standard analysis 
method using a spherical NFW is equivalent to putting 5-function 
priors on both axis ratios at a = & = 1. This is clearly a very strong 
prior, and its application to lensing halos with significant triaxial 



structure was shown in CK07 to lead to errors in mass and concen- 
tration estimates of factors of up to 2, in some cases in the same di- 
rection (i.e. increased mass and concentration), causing potentially 
misleading disagreements with the inverse C — M relationship pre- 
dicted in ACDM simulations. 

The opposite of this strong prior is a very weak flat prior on 
both axis ratios truncated at some small value. Here, as a represen- 
tative of this class of prior, we impose 

0.4 s; b s; 1.0 and 

0.3 s; a/6 s; 1.0, (12) 

giving an effective prior on the minor axis ratio 0.12 ^ a ^ 1.0, 
and refer to this herein as the Flat prior. This is a very general and 
weak prior, but, as with the spherical case, does not reflect our true 
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Figure 7. Af200 ~ C* contours and 9 — 4> contours shaded by the axis ratios and, for the first pairing, orientation angles for a triaxial NFW fitted to weak 
lensing by an extreme prolate halo a=b=0.4 in a Line of Sight orientation, under a Flat prior 



prior beliefs about the axis ratio distribution of triaxial halos. We do 
not use the weakest possible prior, 0.0<a = 6^1.0, because we 
find that allowing unphysically tiny axis ratios makes convergence 
of the MCMC runs very difficult. However, the Flat prior could be 
extended to even lower axis ratios and the length of the MCMC 
runs extended, if the problem at hand demands it. 

An intermediate choic e is to use d i stribut ions from CDM sim- 
ulations, such as those of Ishaw et al. I ( l2006h , as the prior distri- 
bution on axis ratios. We fit polynomial functions to the distribu- 
tions of b and a/b, the intermediate:major axis ratio and the mi- 
nor:intermediate axis ratio - the exact functions used are give in 
Appendix [Cl- and the distributions as used are plotted in Figure[T] 
Note that the distribution for b peaks at about 0.8 and that for a/h 
peaks at about 0.9, giving a population, as found in several other 
simulations as well (see e.g. Bett et al. (2007)), skewed slightly to- 
ward prolate halos. We will call this the Shaw prior herein. 

Additionally, one might wish to place a prior on the mass of 
the lensing halo, based on a simulated or observed mass function. 
For simplicity, we adopt a simple exponential prior 

p(M2oo) oc exp (-A/2oo/1O'''M0) (13) 

to understand the strength of this class of prior compared to others, 
and refer to this as the Mass prior. 

There are many other possible choices of prior, including a 



more general prior on the axis ratios tending towards less extreme 
values than the Flat prior, a prior based on the C — M relationship 
predicted by ACDM, a prior taking into account selection biases 
due to lensing efficiency (some halo shapes and orientations - i.e. 
prolate line of sight - are much better lenses than others), or priors 
based on knowledge about a specific lensing system gained from 
outside data or unusual characteristics of the system. In addition, 
many of these priors may be used in combination. Exploring the 
impact of priors is relatively computationally cheap because they 
can often be applied to existing MCMC chains that have been run 
under more general priors; a post-processing conversion can be per- 
formed as long as everywhere the new prior is non-zero there was 
non-negligible probability under the original prior. Thus, the Shaw 
and Mass priors can usually be applied after the fact to runs under 
the Flat prior. By contrast, the Spherical prior cannot, as there is 
too little relative probability and therefore far too few samples at 
a = b = 1 under the Flat prior. 

The choice of prior is always a balance between imposing 
enough external constraints to gain usefully specific parameter es- 
timates without falsely over-constraining the problem with priors 
too strong or simply wrong for the problem at hand. We explore 
the behaviour of the triaxial NFW fits under several of these priors 
in the following section, and return to this important question in the 
Discussion. 



An MCMC Fitting Method for Triaxial Dark Matter Haloes 9 



o 




.5 



Figure 8. M200 — C contours for a triaxial NFW fitted to weak lensing by 
an extreme prolate halo a=b=0.4 in a line of sight orientation under several 
priors. The black star gives the true parameters of the underlying lens. The 
black contours are the result for a triaxial NFW fit with both axis ratios 
fixed at the the true values a = b = 0.4; as expected, they enclose the true 
parameter values for mass and concentration. 

4 RESULTS AND ANALYSIS 
4.1 Fitting to highly triaxial halos 

To begin to understand both the need for this triaxial fitting method 
and its behaviour under different priors, we examine four cases of 
highly triaxial halos studied in CK07. They are prolate and oblate 
halos of M200 = lO^^Af© and C = 4 with axis ratios (a = 6 = 
0.4) and (a = 0.4; b — 1.0) respectively, with effective spheri- 
cal virial masses and concentrations (m2oo ~ 0.90 x lO^^M©, 
Cph = 3.86) and (maoo = 0-91 x lO^^Mo, Csph = 3.87) 
oriented in Line of Sight and Plane of the Sky orientations; Fig- 
ure |2] shows the convergence contours for these four lensing con- 
figurations. CK07 showed that these four lenses, if fit with spherical 
models using a standard maximum likelihood technique give mean 
parameter values that differ significantly from the true values of the 
lens; 

• Prolate LoS; M200 = 1.34 x 10^^Mq;C = 6.0 

• Prolate Plane; M200 = 0.75 x lO^'^M©; C = 3.4 

• Oblate LoS; M200 = 0.72 x IO^'^Mq; C = 2.7 

• Oblate Plane; M200 = 1-05 x 10^^ Af©; C = 5.0. 

Figure[3]shows the resulting M200 — C contours fitting a triax- 
ial NFW to a single lensing realization through each of these four 
cases applying a Flat, Shaw, Mass, joint Shaw and Mass (Multi), 
and Spherical prior. Figure |4] plots the Id marginalized parameter 
distributions. As expected the broader the prior, the broader the er- 
ror contours. In these extreme cases the Flat prior performs best; 
this is unsurprising in comparison to the Shaw prior, which has 
low probabilities of axis ratios so extreme as those of these halos, 
and certainly in comparison to the spherical prior. The Mass prior 
proves weak, reducing the contours only slightly from the Flat prior 
case, in the expected direction away from high mass solutions. 

The four halos show trends consistent with those expected 
from previous work in their relative masses and concentrations. 
For comparison. Figure [5] shows the effective spherical mass and 



concentration m2oo — Csph contours, and their Id marginalized 
parameter distributions, obtained fitting a triaxial NFW under all 
priors to the prolate and oblate Line of Sight cases. The contours 
and parameter distributions appear very similar to their fully triax- 
ial counterparts, and exhibit the same trends of size and accuracy. 

None of the contours or Id parameter distributions for the Pro- 
late LoS case include the true value of the underlying halo in either 
the triaxial or effective spherical parameterisation. To better under- 
stand this. Figure |6] plots the full posterior probability distribution 
of the triaxial model under a Flat prior, and Figure |7] plots several 
3D plots, colouring the A/200 — C contours and the 6 — (j) contours 
by the axis ratios and, for the first pairing, orientation angles. These 
allow better understanding of what degeneracies exist between the 
parameters, and thus what parts of parameter space are opened by 
allowing extreme axis ratios or particular halo orientations. 

First, we see that the failure of even the very general Flat prior 
contours to enclose the true halo value is due to the intrinsic priors 
on the orientation angles, which make the Line-of-Sight orientation 
of the lensing halo highly unlikely (for this lensing configuration 
9 = Q, which under random orientation has probability approach- 
ing zero). This can be understood intuitively by imagining oneself 
at the centre of this "rugby-ball" halo; there are only two places on 
the sky above where observers can look straight "down the barrel" 
of the halo - those two special points situated directly at the pointed 
ends of the ball. Conversely, there are many places on the sky where 
observers can take in identical side views of the halo ~ all the many 
points radiating out from the circumference of the fattest part of the 
"ball." 

However, despite its statistical unlikelihood, this geometry is 
still important as a limiting case, as it is the strongest lensing con- 
figuration possible for a halo of a given mass and minor axis ra- 
tio. Thus, though halos in such orientations are uncommon, they 
are strongly favoured in lensing-selected samples. Further, they are 
particularly dangerous in that they show little (in the symmetric 
limit shown here - no) ellipticity on the sky, and so are likely to 
be treated as spherical if triaxial modelling is used only selectively. 
In some analyses of very powerful lenses, a more complex prior, 
taking into account lensing efficiency in the prior distribution of 
axis ratios and orientation angles, may be required to capture the 
true posterior distribution. Looking at panels 1-3 in Figure |6] high- 
lights this, as we see that it is low axis-ratio, low 6 (close to Line of 
Sight), low mass solutions that give the models closest to the true 
lens; those that would be favoured by a strong lensing efficiency 
prior Figure[8]shows the resulting contours in the completely unre- 
alistic scenario in which the true underlying axis ratios are known, 
here set to be a = 6 = 0.4, while the mass, concentration, and 
orientation angles remain free. The contours now contain the true 
model values for mass and concentration, confirming that the tri- 
axial fitting routine behaves predictably and correctly in limits of 
both maximum and minimum knowledge about the underlying lens 
geometry. 

Figures [9] and [TO] show the posterior probability distribution 
for the triaxial model under the Shaw and Mass priors, omitting the 
orientation angles for brevity. As expected, the axis ratio distribu- 
tion is significantly constrained under the Shaw prior; this leads to 
tighter Af2oo ^ C contours that overall favour lower masses and 
concentrations, but are completely contained within those obtained 
under the Flat prior This indicates that some very low mass, low 
concentration models have also been lost. Clearly, the Shaw prior, 
which has very little probability in the region where this halo ac- 
tually sits in parameter space, is not a good prior to impose in this 
case. The Mass prior, while weak, does move the contours in the di- 
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rection of the true model, eliminating some higher mass cases while 
keeping the very low mass, low concentration cases. This suggests 
that in cases in which there are reasons to suspect a particularly 
advantageous lensing geometry - perhaps because a studied lens is 
one of the strongest in our observable universe, or dynamical stud- 
ies suggest elongation along the line of sight - a prior favouring 
lower masses may be an elegant and clear way to account for that 
prior belief, in place of a direct lensing efficiency prior. 

When applying this MCMC triaxial NFW fitting method (or 
any other method) to real lensing data from galaxy clusters, there 
may be some uncertainty in the exact location of the cluster cen- 
tre. This may be accounted for if necessary by including the centre 
position of the halo model as an unknown in the fit, increasing the 
number of free parameters by two for all models. However, we find 
that small errors in the fixed position of the cluster centre lead to 
only small changes in the most probable parameters recovered by 
our MCMC method, and that those changes are of the same scale 
under all priors, including at the extremes the highly triaxial Flat 
prior and the over-constrained Spherical prior The errors induced 
by fixing the cluster centre are therefore no more a problem in a 
triaxial analysis than in any other, and we neglect them herein in 
our assessment of the behaviour of our new triaxial NFW fitting 
method. 

4.2 Fitting to a population of Iialos 

We now move from looking at the behaviour of the triaxial NFW 
fitting technique on individual, highly triaxial halos, to its be- 
haviour across a physically-motivated population of triaxial halos. 
We use a standard Monte Carlo technique to choose 100 triaxial 
NFW halos with M200 = lO^ Mp and C = 4 and axis ratios 
drawn from the distributions in IShaw et al. I ( I2OO6I) . The parame- 
ter distributions of the population are shown in Appendix |Bj the 
mean effective spherical parameters for the population are similar 
to the constant triaxial values: < m2oo >= 0.99 x 10^^ A/0 and 
< Csph >= 3.98. We randomly orient the hundred halos, lens 
through them as described in i]2.3l and carry out MCMC fits of the 
triaxial NFW model under the various priors previously described. 
FigurefTTIshows the resulting most-probable triaxial parameter dis- 
tributions for the population under a Flat prior and compares them 
with a Spherical prior typically employed in lensing analyses. Fig- 
ure [12] shows the same results under a Shaw prior, which exactly 
matches the true shape distribution of the population. 

Under the Flat prior, the distribution of mean best-fitting pa- 
rameters is skewed generally towards higher values of concentra- 
tion and mass compared to that under the Spherical prior. Despite 
the underconstrained nature of the problem, the recovered axis ra- 
tios show a significant dependence on the lensing data, as uniform 
priors for the axis ratio distributions become distributions peaked 
towards higher values of a and b. Under the Shaw prior, the triaxial 
models give mean best parameter values very similar to those under 
the Spherical prior. This reflects the similar shapes of the posterior 
probability confidence contours (e.g. the similar 'banana'-shaped 
M — C contours) under these two priors. 

The upper limits of the intermediate and minor:intermediate 
axis ratios are strongly constrained by the visible projected elliptic- 
ity on the sky (b or a/b must be less than or equal to q), while a soft 
lower limit is imposed by the fact that there are only a few geome- 
tries that allow for axis ratios significantly less than the observed 
ellipticity to be hidden in projection. 

Note that throughout we quote mean best-fitting parameter 
values as opposed to peak probability values; this gives parame- 



Table 1. The percentage of triaxial NFW fits to a population of 100 lens- 
ing halos (with A/200 = 10^^ M©, C = 4 and axis ratios drawn from 
the ACDM-motivated Shaw axis ratio distribution) for which the true mass 
and concentration values fall within the 68% and 95% confidence contours, 
under various priors on the halo parameters. 



Prior 


68% 


95% 


Flat 


86 


99 


Shaw 


66 


94 


Axis 


70 


96 


Mass 


86 


99 


Spherical 


53 


81 




Effective Spherical Parameterisation 




Flat 


84 


99 


Shaw 


61 


89 


Spherical 


56 


82 



ter values further away from the extremes of the distribution. This 
is seen, for example, looking back at the ID parameter distribu- 
tions for the extreme halos in Figured although the probability 
distribution for b peaks at 6 = 1 for both the Prolate and Oblate 
LoS cases, the mean best-fitting values are b = 0.89 and b = 0.88 
respectively. 

Figure [T3] shows A/200 — C confidence contours for a typical 
triaxial halo from the population, with axis ratios {a = 0.76, b = 
0.85}. Note that opening up the axis parameter space generally 
leads to many models with high mass and high concentration com- 
pared to the over-constrained spherical model, and only a few lower 
mass, lower concentration models. This is the general behaviour of 
the posterior probability distribution for all halos, and indicates that 
neglected triaxiality in lens models cannot explain a whole popula- 
tion of massive halos with high concentrations within the ACDM 
paradigm. It is only in a small fraction of orientations that highly 
triaxial halos lead to highly efficient lensing and thus to overesti- 
mates of mass and concentration under a spherical prior This is 
seen in the second panel of Figure [T3] plotting the A/200 — C con- 
tours under the Flat prior coloured by minor axis ratio a. Most of 
the small-a models are in the upper right corner of the distribu- 
tion, but there is also a narrow band of them forming the lower left 
boundary of the posterior distribution, representing the highly ef- 
ficient triaxial lens orientations that lead to overestimates of mass 
and concentration under standard spherical analysis methods. Fur- 
ther note that the familiar degeneracy between mass and concen- 
tration is still evident for each range of minor axis ratio, with the 
allowance of more extreme axis ratios adding additional parallel 
lines of likely models, as indicated by the colour banding. 

The top half of table □ gives the number of cases in which 
the true triaxial model value falls within the 68% and 95% confi- 
dence contours under each prior. The bottom half of the table gives 
the equivalent results for the Flat, Shaw, and Spherical prior using 
the effective spherical parameterisation for both the fitted triaxial 
models and the triaxial lenses they are fit to. Most strikingly, un- 
der the triaxial parameterisation the numbers correspond very well 
to the predicted statistics when the correct prior is employed; un- 
der the Shaw prior indeed 66 of the halos include the true value in 
their 68% contour, and 94 include the true value in their 95% con- 
tour! This is very good evidence that the MCMC fitting technique 
behaves statistically as expected, in that the errors under a correct 
prior are the true errors due the full uncertainties in the problem: 
the statistical uncertainty in the shapes of the background galaxies, 
and the geometric uncertainty regarding the line-of-sight structure 
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Figure 9. The full posterior probability distribution for a triaxial NFW fitted to weak lensing by an extreme prolate halo a=b=0.4 in a Line of Sight orientation, 
under a Shaw prior. The black contours give the 68% and 95% confidence contours of the marginalized posterior probability distribution, and the shading 
shows the marginalized likelihood values, smoothed over 20 bins (rather than 40 as in other figures) for increased clarity. The solid (dotted) lines in the 
unshaded distributions at the end of each row give the 1-parameter marginalized probability (likelihood) distributions. 



of the lensing halo. As expected, the very general Flat prior gives 
larger contours that in a larger number of cases contain the true val- 
ues. The only slightly stronger Mass prior behaves similarly, and 
the heavily over-constrained Spherical prior gives small contours 
that far too often exclude the true parameter values. 

Since we will never know the exact axis ratio distribution of 
the galaxy cluster population, we test the response of the method 
applied to the same Shaw population of halos usi ng a slightly differ - 
ent physically motivated prior taken loosely from lBett et al.l ( [2007h , 
plotted in Figure [T] and called herein the Axis prior. We define 
it simply by two normalised Gaussians, with p{b) having mean 
^ = 0.8 and standard deviation a — 0.125 and p{a/b) slightly 
narrower with /i = 0.85 and a = 0.1. The result is encouraging; it 
is almost as good an estimator of both the mean parameter values 
across the population and the error contours on those parameters. 
Thus, a slight mismatch between physically motivated priors and 
the true axis ratio distribution will not significantly decrease the 
accuracy and usefulness of this method. 

Using the effective spherical parameterisation, fits employ- 
ing the triaxial Shaw prior still statistically outperform the over- 
constrained Spherical prior, demonstrating that the triaxial model 
with a well-chosen prior better recovers even effective spherical pa- 



rameters than does an over-simple spherical model. However, fits 
under the Shaw and Axis priors using the full triaxial parameter- 
isation are the overall best statistical performers, supporting our 
choice to use the fully triaxial parameterisation whenever possible. 

The top half of Table|2]gives the mean triaxial mass and con- 
centration values recovered across the population (fiducial values 
Af200 = 1-0 X IO^^A/q and C = 4); the bottom half gives 
the means of their effective spherical counterparts under the Flat, 
Shaw, and Spherical priors (fiducial mean values m2oo ~ 0.986 x 
1O^^M0 and Csph = 3.98). Interestingly, and in line with earlier 
findings in CK07, the spherical model serves, across this realistic 
triaxial halo population, as a very good estimator of the mean triax- 
ial and effective spherical mass and concentration. While the fact 
that so often the contours under the Spherical prior do not contain 
the true parameter values makes it clear that the model is a bad 
choice for individual halos (it does not accurately constrain mass 
and concentration simultaneously, and the error estimates it pro- 
vides are signficantly too small), this result suggests that if only 
mean values are needed, the spherical model may be adequate. 

However, in the context of measuring a mass function or un- 
der poorer observing conditions where the error contours are more 
asymmetric, its inability to constrain individual halo models or their 
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Figure 10. The full posterior probability distribution for a triaxial NFW fitted to weak lensing by an extreme prolate halo a=b=0.4 in a Line of Sight orientation, 
under a Mass prior. The black contours give the 68% and 95% confidence contours of the marginalized posterior probability distribution, and the shading 
shows the marginalized likelihood values, smoothed over 20 bins (rather than 40 as in other figures) for increased clarity. The solid (dotted) lines in the 
unshaded distributions at the end of each row give the 1-parameter marginalized probability (likelihood) distributions. 



errors accurately make the spherical model most likely inadequate. 
This question will be further investigated in an upcoming paper 
by the authors. As expected, the triaxial fitting method under the 
Shaw prior is also a good estimator of the mean population values 
in both the triaxial and effective spherical parameterisations, and 
crucially also estimates the errors accurately. The Flat prior and the 
closely related Mass prior result in overestimates of mass and con- 
centration, as expected given the shape of the posterior probability 
distribution (many high-M, high-C models). 



5 DISCUSSION & CONCLUSIONS 

Some may argue that it is foolhardy to fit a density profile model 
that fundamentally cannot be fully constrained with available data. 
Indeed, in many cases calculating the Bayesian evidence would 
likely favour a Spherical over a Flat prior, due to a significant de- 
crease in available parameter space coupled with a relatively small 
decrease in likelihood values under the Spherical prior. However, 
unlike in cases of fitting multiple halos to account for substructure, 
or in another context, adding additional parameters to cosmolog- 
ical models, we know a priori from physical observations that a 
triaxial model is a better model than a spherical model - we see 



Table 2. The mean most-probable parameter values resulting from fitting 
a triaxial NFW to a population of 100 lensing halos (with M200 = 10^^ 
Mq , C = 4 and axis ratios drawn from the ACDM-motivated Shaw axis 
ratio distribution), under various priors on the halo axis ratios and mass. 



Prior 


Af2Oo[lOi5M0] 


C 


Original Population 


1.00 


4.0 


Flat 


1.12 


4.7 


Shaw 


1.04 


4.2 


Axis 


1.04 


4.3 


Mass 


1.07 


4.7 


Spherical 


1.02 


4.1 


Effective Spherical Parameterisation 


"1200 


C sph 


Original Population 


0.986 


3.98 


Flat 


1.06 


4.6 


Shaw 


1.03 


4.2 


Spherical 


1.02 


4.1 



clearly in non-parametric lensing mass maps that galaxy clusters 
have significant levels of ellipticity and are not spherical. The pa- 
rameter estimates and uncertainties resulting from spherical mod- 
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Figure 11. Parameter distributions resulting from fitting the triaxial NFW to a randomly oriented population of 100 tensing halos with M200 = 10 Mq, 
C = 4 and axis ratios drawn from the Shaw prior. The results under the very general Flat prior are compared to those under a Spherical prior The top left 
panel plots the mean best-fitting concentration and mass values for each lens in the population under the two priors; the top light panel plots the returned axis 
ratio distributions (histograms) under the Flat prior compared to the prior distributions themselves (p(b) solid line, p(a/b) dashed line, p(a) dotted line); the 
bottom panels plot the distributions of mass and concentration obtained under each prior 



els, or 2D elliptical models, that assume a (5-function prior on at 
least one axis ratio, therefore intrinsically claim a level of certainty 
we know we do not possess. Making our best choice of prior on 
the axis ratios and fitting triaxial models is the more physically true 
approach, and gives the best estimate of parameter values and er- 
rors that can be obtained from tensing data alone. The large size of 
these contours emphasises the need to combine tensing data with 
other data sources, be they SZ, dynamical, or X-ray, because the 
dependence of all tensing effects - strong, weak, intermediate - on 
projected mass fundamentally constrains the precision and accu- 
racy with which galaxy cluster masses and concentrations can be 
obtained. 

Further, there is no way to know for a single tensing cluster 
how important triaxiality is without first carrying out a full triax- 
ial analysis and comparing its results to those from a spherical fit; 
Figure [14] helps illustrate why this is the case. The four panels plot 
the distributions of two observable properties of the triaxial halos 
of the tensing population of Section |42l for which the 68% and 
95% confidence contours obtained fitting an NFW under a Spheri- 
cal prior do not contain the true halo values. These observables are 
the axis ratio of the projected isodensity contours q of the lens and 
the maximum likelihood value of the NFW fit under a Spherical 
prior. It might be hoped that either a high visible ellipticity (low q) 
or a low likelihood value would indicate that the spherical model is 
a bad fit; however, as seen in the figure, the distributions show no 
bias compared to the distribution of the entire tensing population, 
and thus can provide no indication of the adequacy of the spheri- 
cal fit. This can be understood by example by referring back to the 




Figure 13. The left panel shows the Af200 — C 68% and 95% confidence 
contours, under various priors, for a triaxial NFW fit to a triaxial halo typical 
of those found in simulations, with Af200 = 10^^ and C = 4, {a = 
0.76, b = 0.85}. The right panel plots the M2Q0 — C distribution shaded 
by the minor axis ratio a for the same halo under a Flat prior. 

extreme halo shapes and orientations of Section HTl the highly tri- 
axial prolate halo in a Line of Sight orientation will have both a low 
ellipticity {q ^ 1) and a very high maximum likelihood value (it is 
an efficient lens with circular projected isodensity contours and so 
will be very well fit by a spherical model), but the most-probable 
parameter values and their errors obtained from fitting a spherical 
NFW are very inaccurate. Thus, it is necessary to fit a triaxial NFW 
using the best prior available to every lens, as there is no way to de- 
termine beforehand how triaxial any given halo might be. 

A great advantage of this MCMC method is the ease with 
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Figure 12. Parameter distributions resulting from fitting the triaxial NFW to a randomly oriented population of 100 lensing halos with M200 = 10 Mq, 
C = 4 and axis ratios drawn from the Shaw prior. The results under the Shaw prior (which matches the distribution of the lens population) are compared to 
those under a Spherical prior. The top left panel plots the mean best-fitting concentration and mass values for each lens in the population under the two priors; 
the top right panel plots the returned axis ratio distributions (histograms) under the Shaw prior compared to the prior distributions themselves (p(b) solid Hne, 
p(a/b) dashed line, p(a) dotted line); the bottom panels plot the distributions of mass and concentration obtained under each prior. 



which other data sets and constraints may be incorporated. Whether 
through extra terms in the likelihood function (though such inclu- 
sions require a careful analysis, not seen yet in the literature for the 
combination of many data types, of the proper weights to give data 
types with very different scales and sources of errors), or through a 
prior (e.g. a prior on the cluster mass taken from the best estimate 
and errors of a dynamical study). 

The future of understanding cluster structure lies in the com- 
bination of various data types, as all available methods are limited 
either by fundamental constraints, such as lensing, or by the avail- 
ability of data, as in dynamics, in which highly simplified dynami- 
cal models are required by the relatively small numbers of position- 
velocity pairs. Methods such as this, and those currently being de- 
veloped by e.g., IJulloetalJ ( l2007h andlFeroz & Hobsod jlOOSh . are 
a crucial element of the next generation analysis toolbox for galaxy 
cluster studies. 



5.1 Selecting better priors 

Placing well-motivated priors on the shapes of galaxy clusters is 
crucial to future work modelling the most massive structures in the 
universe, and understanding both their characteristics as individu- 
als and a population. While structure formation simulations pro- 
vide a good starting point for such priors, using their predictions 
does of course bias all results towards agreement with those simu- 
lations. One good alternative is to use a distribution of 2d axis ra- 
tios observed in a sample of galaxy clusters, preferably from mass- 
sensitive methods (e.g. lensing mass reconstructions) and from this 



construct a 3d shape prior. For this, even the simplest elliptical lens- 
ing models would be adequate, as CK07 showed that even a Sin- 
gular Isothermal Ellipsoid model consistently recovered projected 
axis ratios accurately, even if the true lensing profile was NFW (this 
is not true for other model parameters). 

While the necessity and importance of these priors may be 
discomforting, using standard techniques and simpler models is 
simply disguising the problem, as such models contain highly- 
restrictive hidden priors. Bayesian techniques such as the one pre- 
sented in this paper are simply tools that allow us to better un- 
derstand the true constraints we can place on physical models, not 
solutions in themselves to physical problems. The broad posterior 
probability distributions for the triaxial NFW profiles fit to simu- 
lated lensing data in this paper indicate the weakness of our current 
constraints, and emphasise the need for focus on the careful com- 
bination of complementary data types to further constrain galaxy 
cluster structure models. To this end, this Bayesian MCMC triaxial 
NFW fitting method provides, through the prior probability func- 
tions, a statistically robust and straightforward way to combine con- 
straints from data types with very different error properties. 

5.2 Application 

This MCMC triaxial NFW fitting method, using fully tested look- 
up tables to significantly speed up the calculation of the triaxial 
NFW lensing quantities, can be implemented for a standard weak 
lensing data set on a single fast processor in about 1 day, making 
it fully feasible for use across even the largest existing lens sur- 
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• allows great flexibility in the inclusion of constraints from dif- 
ferent observations and analyses, including strong lensing, dynam- 
ical studies, and S-Z and X-ray derived mass models; 

• demonstrates that while triaxiality can explain rare cases of 
overconcentrated galaxy cluster halos within the ACDM paradigm, 
it cannot explain a population-wide trend of over-concentration; 

• should be widely applied in weak lensing analyses to better 
determine cluster parameters and errors that represent the true (and 
limited) extent of our knowledge of the distribution of matter in 
galaxy cluster lenses. 




Scaled max lil^elihood Scaled max lil^elitiood 



Figure 14. The right (left) top and bottom panels show, respectively, the dis- 
tributions of the projected axis ratios q and the maximum likelihood values 
obtained from fitting a spherical NFW, for the halos of the triaxial lensing 
population of Section 1431 for which the 68% (95%) confidence contours 
recovered fitting an NFW under a Spherical prior do not enclose the true 
triaxial halo values. 

veys. To begin this process, in a companion paper Corless, King, & 
Clowe (in preparation), we apply this method employing a range of 
priors to galaxy clusters Abell 1689, Abell 1835 and Abell 2204. 

In applying this method to observational data it will be im- 
portant to account for associated and line-of-sight structures in the 
lensing field. Galaxy clusters form at the intersections of filaments, 
so it is not surprising that spectro scopic studies ofte n reveal corre- 
lated structures in their fields(e.g. lLokas et al.l (2006)) that can bias 
our es timation of cluster parameters if neglected (King & CorlessI 
1 I2OO7I) ). Unaccounted for large scale structure along t he line of 
sight c an also have a severe impact on lensing analyses jHoekstrj 
l l2003l) ).leading to a factor of ^ 2 incr ease in errors on cluster mass. 
As demonstrated by 'Po delsonl ( l2004j ). this can be somewhat miti- 
gated by including a noise term for large scale structure in the lens- 
ing analysis. Building on the greater coherence scale of large scale 
structure noise compared with the noise associated with the intrin- 
sic ellipticity dispersion of galaxies, Dodelson notes that the errors 
on cluster mass can be reduced by ^ 50% for wide-field data. 

This MCMC triaxial fitting method will be very useful in 
the analysis of current and future surveys such as LoCuSS, PAN- 
STARRS and DES. 

5.3 Summary 

We here present an MCMC method for fitting triaxial NFW models 
to weak lensing data that 

• consistently and accurately returns parameter and error esti- 
mates representing the true uncertainties of the problem; 

• includes Bayesian priors on the halo geometry, chosen from 
simulations or observations; 
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potential $ (commas indicate differentiation): 



71 = 

72 = 



2 ($,xx + '^,yy) 



(AlO) 
(All) 
(A 12) 



These derivatives are calculated as functions of integr als of the 
spherical convergence Hi{C,) (see e.g. iBartelmannI ( Il996h for a full 
treatment of weak lensin g by a spherical NFW profile) following 
the method of lSchramml i 1990) and Keetori c2001i) . normalised by 
a factor of l/\ff from Equation lA4l (see OLS for the derivation of 
this normalisation) 



= 2qX^Ko + qJo, 



2qY'K2 4 
2qXYKi, 



qJ\ 



(A 13) 
(A14) 
(A15) 



APPENDIX A: LENSING THROUGH TRIAXL4L HALOS 

Following OLS, the triaxial halo is projected onto the plane of the 
sky to find its projected elliptical isodensity contours as a function 
of the halo's axis ratios and orientation angles {6, (j)) with respect 
to the the observer's line-of-sis htQ The elliptical radius is given by 



.2 

C = — + — 

9| 4 



(Al) 



where (X, Y) are physical coordinates on the sky with respect to 
the centre of the halo, 



2 



V 



where 



■2 i <? (? 
f — sin 6^ — cos + -T- sin ( 



+ cos 



(A2) 
(A3) 

(A4) 



and 

A 



2 a / . 2 , 2 

cos 6' I ^ sm (p + jTT cos ( 



B — cos S sin 2(7!> ^ — ttt 

b^ 

2 2 
n ■ 2 , C 2 , 

C = — sm H -cos 0. 

0^ 



The axis ratio q of the elliptical contours is then given by 

and their orientation angle on the sky by 
X > gy). 



2 2 

+ ^l2sin'e, (A5) 



(A6) 
(A7) 

(A8) 
(A9) 



1 _i B 
* = - tan - 

2 A~C 

Here we diverge slightly from OLS's treatment as we are in- 
terested not in deflection angles but in the lensing shear and con- 
vergence, both combinations of second derivatives of the lensing 



where 

A"u(x,y) 
Jn(x,y) 

and 

CM' 



VTJo [l-(l-g2)u]"+i/2 



77 7o [i-(i-g^H"+^/^ 



Y' 



du. 



du. 



qx \ 1 — (1 — q^)u 



(A16) 
(A17) 

(A 18) 



Note that our radial variable ^ appears different from Keeton's 
^ because it is defined in terms of two axis ratios qx and qy rather 
than one q: ( — £,/qx- This reflects a dependence on the 3D struc- 
ture of the cluster; for example, extended structure along the line of 
sight decreases qx and thus increases the convergence and shear at 
a given {X, Y). 



APPENDIX B: PARAMETERS OF THE LENSING 
POPULATION 

The 100 halos that make up the lensing population of Section [42] 
each have triaxial virial mass Af200 = 1.0 X 10^^ Mq and con- 
centration C — 4.0. Their axis ratios are drawn fro m the distri- 
bution s found in the structure formation simulations of IShaw et al.l 
( l2006l) , and their effective spherical virial masses m2oo and con- 
centrations Csph are calculated numerically as a function of their 
triaxial shape. The distributions of axis ratios a and b and the ef- 
fective spherical parameters m2oo and Csph are plotted in Figure 



^ Although we set c = 1, we keep c as a variable in our notation for 
consistency with OLS. 



An MCMC Fitting Method for Triaxial Dark Matter Haloes 1 7 




Figure Al. The distribution of axis ratios a and b and effective splierical 
mass m200 and concentration C^ph for 'he lensing population of Section 

m 
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APPENDIX C: PRIORS FROM SIMULATIONS: SHAW PRIOR 

We define the Shaw prior, plotted in Figure[T] by fitting polynomials to the data points of Figure 14 in lShaw et alj j2006h : 
if6<0.5 
[l.6329fe^ - 7.97756* + 9.3414fo^ - 6.6558fe^ + 2.29646 - .3088] x 10^ if 0.5 6 < 1.0 



p{h) 
and 



'(f) 



if f < 0.65 



5.76647 (f) - 2.459265 (f) + 42.3154 (f) - 37.2765 (f) if 0.65 f ^ 1.0. 



17.4650 (f )^ - 4.00238 (f ) + .32462 



X 10* 



This paper has been typeset from a T^X/ ETbX file prepared by the author. 



